Transport through two interacting resonant levels connected by a Fermi sea 
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We study transport at finite bias, i.e. beyond the linear regime, through two interacting resonant 
levels connected by a Fermi sea, by means of time-dependent density matrix renormalization group. 
We first consider methodological issues, like the protocol that leads to a current-currying state and 
the characterization of the steady state. At finite sizes both the current and the occupations of the 
interacting levels oscillate as a function of time. We determine the amplitude and period of such 
oscillations as a function of bias and extension of the Fermi sea. In particular, the occupations on the 
two dots oscillate with a relative phase which depends on the distance between the impurities and on 
the Fermi momentum of the Fermi sea, as expected for RKKY interactions. Also the approximant 
to the steady-state current displays oscillations as a function of the distance between the impurities. 
Such a behavior can be explained by resonances in the free case. We discuss finally the incidence of 
interaction on such a behavior. 

PACS numbers: 72.10.Fk, 73.23.-b,73.63.Kv, 75.40.Mg 



I. INTRODUCTION 

The study of quantum transport across nanostructures 
has been the subject of intense theoretical and experi- 
mental attention for decades. One of the most intensively 
studied systems is that of quantum dots, both because 
of their great experimental versatility and because they 
unveil an extremely rich physics, as exemplified by the 
Kondo effect^ in quantum dots 2 ^. When considering a 
system of two quantum dots, a further interesting phe- 
nomenon emerges, the Rudcrmann-Kittcl-Kasuya-Yosida 
(RKKY) interaction^. It describes the indirect interac- 
tion between two magnetic impurities mediated by the 
electrons of the surrounding Fermi sea, and is charac- 
terized by oscillations related to the Fermi wavevector. 
The competition between the RKKY interaction and the 
Kondo effect was studied in the frame of numerical renor- 
malization group£, and conformal field theory^ 7 -. An ex- 
perimental realization with two quantum dots coupled by 
a Fermi sea was meanwhile reported^. 

Recently, a great deal of progress was achieved towards 
the theoretical determination of steady-state transport 
properties focusing on a quantum dot described by the 
interacting resonant level model (IRLM)2r— , that con- 
sists of spinless fermions with a nearest-neighbor re- 
pulsive interaction for the sites adjacent to the dot. 
This model was studied with several theoretical tech- 
niques, ranging from integrable field theories and Bethe 
Ansatz (see Boulat et ali£ and references therein), func- 
tional renormalization groupiSr— , real-time renormaliza- 
tion group^, to density-matrix renormalization group 
(DMRG) techniques2*i^i2 . These works provide the I-V 
characteristics out of equilibrium at finite bias and up 
to large values of the interaction^, and a detailed knowl- 
edge of the relaxation dynamics ^ 12 ' 14 in the regime of 
small interaction, including also the incidence of finite 
temperatures^. The shot noise and the full counting 
statistics have been studied by means of exact diagonal- 
ization-^ (in the free case) , DMRG and thermodynamical 



Bethe Ansata!^. Such an attention on a model that 
arguably cannot be experimentally realized in an elec- 
tronic system is due to the fact that, in contrast to the 
Anderson impurity model, the important energy scales of 
the problem are accessible and controllable in numerical 
simulations, avoiding to deal with the Kondo scale, that 
requires high resolution in energy. 

In contrast to the great attention devoted to the one 
impurity case, little is known about the case with more 
impurities2i~— . In particular, to the best of our knowl- 
edge, the case of two IRLs separated by a Fermi sea under 
a finite bias awaits still a theoretical treatment. Here we 
consider two leads modeled as tight-binding chains with 
uniform hopping, coupled to two quantum dots inter- 
acting with their nearest-neighboring sites and a Fermi 
sea in between, focusing on the dynamics of the system 
when it is taken out of equilibrium with the application 
of a finite bias. The set-up is that of a quantum quench, 
where the initial state corresponds to the ground state of 
a Hamiltonian, but the time evolution is governed by a 
different (time independent) one. We considered two dif- 
ferent protocols, where the bias is included either in the 
initial or in the final Hamiltonian. We discuss also the 
characterization of the steady-state and the incidence of 
finite-size effects. 

We performed our studies by means of a time- 
dependent DMRG (t-DMRG) simulation 2 ^. This 
method allows to study the time evolution of the sys- 
tem up to intermediate times (~ 4QH/ta, where to is the 
nearest-neighbor hopping between the sites of the leads) 
in a nonperturbative way. The time evolution of the cur- 
rent on each link of the chain and of the particle-density 
on the dots exhibits oscillations that depend not only on 
the bias, as in the single dot case, but also on the exten- 
sion of the Fermi sea. Furthermore the occupations on 
the two dots oscillate with a relative phase which depends 
on the distance between the impurities and on the Fermi 
momentum of the Fermi sea. This can be explained in 
terms of the RKKY interaction. The currents in the sites 
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connecting the quantum dots to the leads show also oscil- 
lations but with a phase shift with respect to the density. 
For the approximant to the steady-state current we find 
that it oscillates as a function of the distance between the 
impurities. In the free case the behavior of the current 
can be understood in terms of resonances that appear 
in the transmission coefficient of a single particle prop- 
agating through the system. We finally show the effect 
of interaction. While it suppresses the resonances found 
in the free case, for strong interaction we find that large 
oscillations of the current as a function of the interimpu- 
rity distance arise, with a periodicity which matches the 
RKKY prediction. 

This paper is organized as follows. Section|H]is devoted 
to the discussion of methodological issues. In particular 
in Sec. Ill Al we define the model, the observablcs and 
the numerical technique. We show the effect of different 
quench schemes and motivate our choice in Sec. Ill B[ In 
Sec. Ill CI we detail how the approximant of the steady- 
state current is obtained and benchmark our DMRG re- 
sults for the one impurity system with those of Boulat 
et al£. Section IIIII displays our results. In Sec. IIII Al 
the time evolution of the occupations and the currents 
is shown and its relation with RKKY interaction is dis- 
cussed. In Sec. IIII Bl we concentrate on the approximant 
to the steady-state values of the current as a function of 
the distance. We consider first the free case, for which we 
establish a connection with the problem of transmission 
of a single particle propagating in the system, and then 
move to the interacting case. In Sec. IIVI we summarize 
our results. 



II. MODELS AND METHODS 



Furthermore, the Fermi sea is described by the Hamilto- 
nian Hp, 



£ (c]c J+ i+H.c.) 



(4) 



In what follows we call the sites at positions Ci = d\ — 1 
and C2 = d-2 + I contacts. The total number of sites of 
the system is given by L, which we take even. If R is 
odd, we choose the position of the dots such that the left 
and the right leads have the same number of sites. If R 
is even the position of the dots is given by (L — R)/2+ 1 
and (L + R)/2 + 1, implying that the left lead has one 
more site with respect to the right one. In Eqs. @ - (HJ) 
we have hj = CjCj, where arc creation (annihila- 

tion) operators for spinless fermions, Uc is the interaction 
coupling the dots and their nearest neighbors, tc is the 
hopping between the dot and its nearest-neighbors. The 
hopping elements in the leads and in the Fermi sea are 
all set to to. Energies are measured in units of to and 
time in units of ft/to- The number of particles in the sys- 
tem is N and we define the average density of particles 
as p = N/L. When not explicitly specified, we assume 
the system at half- filling. We also define L c = R + 1, N c 
and p c = N c /L c as the number of sites, the number of 
particles and the density in the central region (the Fermi 
sea plus the dots), respectively. The system is depicted 
in Fig. [Q 
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A. Hamiltonian and observables 

We study a system characterized by the presence of 
two quantum dots at positions d\ and di separated by a 
distance R = di — d\. The region inbetween harbours a 
Fermi sea. The Hamiltonian of the whole system is given 
by 

-Hchain = #D + Ht + Hp , (I) 

where 

Hb = -ic^-i^i + 4 2 g rf 2 +i + H - c -) 

-*c (4^1+1 + 4 2 -i^2 + H - c -) 



-Uc 



a=d\ ,d 2 r=±l 
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corresponds to the dots and their nearest-neighbors, 
where the interaction is present. The leads connecting 
to the quantum dot are described by the tight-binding 
Hamiltonian H^, 

di-2 L-l 

H T = -to ^ CjCj+i - to c]c 0+ i + H.c. . (3) 

3 = 1 j=d 2 + l 



Figure 1. (Color online) Picture of the model Eq. [JJ for a 
system of L — 14 sites and R — 5. The shaded light blue 
areas indicate the presence of the bias (Eq. [TJl . 



As it will be discussed in more detail in Sec. Ill BL 
we will follow the transport process in the frame of a 
quantum quench, where a given initial state l^o) evolves 
in time under the action of a given Hamiltonian, such 
that the state of the system at a time r is \^(t)) = 
exp(— iHr^^o}- Accordingly, the time-dependent occu- 
pations on each site are given by 



nj (r) = <* (t^I^t)) 



(5) 



The current on each bond connecting nearest-neighbor 
sites can be obtained as: 

h = ^rf J (*(T)\(c]c J+1 - at +1 c,-)l*W> . (6) 

where e is the electron charge, tj is the hopping on the 
bond connecting sites j and j + 1. 

The results presented in this work are obtained with 
the t-DMRG^— . We typically simulate systems with 
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L ~ 100 sites. In order to implement the time evolution, 
we use the Trotter decompositio n 24 ' 25 ' 28 . Our code is 
adaptiv e; 24 ' 25 ' 28 , meaning that the number of states used 
at each time step changes dynamically keeping the dis- 
carded weight below a given threshold. The maximum 
number of states used in our computation is m ~ 1000 
and the discarded weight e is kept below <~ 10~ 7 . In the 
absence of interactions we employ also exact diagonaliza- 
tion. 



B. Quench schemes 

In order to initiate transport processes in the system 
described by Eq. ([T]), a bias AV has to be applied on the 
left and the right lead. It is described by: 



Hr: 




(7) 



As previously discussed for a single impurity^ 6 -, we can 
start with the ground state of -ffchain an d follow the evo- 
lution of the system dictated by a Hamiltonian -ff c hain + 
Hb- Wc denote such a procedure scheme (A). In such 
a scheme, however, the bounded nature of the spectrum 
of a lattice model becomes evident whenever the bias ex- 
ceeds the bandwidth. In that case, there are no states 
available for transport through the system, as shown in 
Fig. [2] (the determination of the currents depicted will 
be discussed in detail in Sec. Ill CI) . It was suggested 
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Figure 2. (Color online) I-V characteristics for a system with 
one (a) and two impurities (b). Black empty circles and red 
full squares refer to quench schemes (A) and (B) respectively. 
All the curves are obtained with L — 100, except for R = 7, 
scheme (A), for which L — 300 sites are used. Data for tc — 
0.8to and Uc = 0.0. The current is shown in absolute value. 



previously 9 -^, that in order to avoid such an artifact of a 
lattice model, the opposite scheme can be used, namely, 
the initial state is the ground state of -f/chain + Hb , and 
the evolution is studied switching off Hb- As shown in 
Fig. O such a quench scheme leads to a saturation of the 
attained current, with similar behavior for a single im- 
purity or two of them with a Fermi sea inbetween. The 
current in scheme (B) saturates at large values of the 
bias, because of the finite bandwidth of the syste m 16 ' 29 . 



For AV smaller than half the bandwidth, both schemes 
lead to the same result. Moreover, for the whole range of 
biases studied in the one-impurity case 9 -, the I-V curves 
can be brought in this way to coincide with analytical 
results from conformal field theory, that correspond to 
the infinite band- width limit. 

In scheme (B) the initial state is characterized by a 
particle imbalance between the left and right lead, due 
to the presence of the bias, and the distribution of par- 
ticles in the central region is not uniform. However, we 
find p c = p if the system is at half filling. In the other 
cases there is a discrepancy which can be controlled by 
performing a finite-size scaling. 

In the rest of the work wc choose quench scheme (B) 
because it avoids the artifact introduced by a bounded 
spectrum. 



C. Time averages 

As already discussed in the Refs. [HI and [13 in the case 
of a single quantum dot, the time evolution of a current 
in a finite system is affected in various ways. On the 
one hand, right after switching the bias on (or off), there 
is a transient time, where the current grows from zero 
to a quasi-stationary state. On the other hand, at long 
times, the current bounces back at the ends of the system. 
In the intermediate quasi-stationary state, periodic vari- 
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Figure 3. (Color online) Finite size scaling of the oscilla- 
tion amplitudes Ij from cosine fits as discussed in main text, 
extracted from the left-contact current I c \. Data refer to a 
system with tc = 0.8to, Uc — 0, 7? = (empty symbols) and 
R = 7 (full symbols). 



ations previously denoted Josephson oscillations 1 ^, due 
to their similarity with the ones in a Josephson junc- 
tion, appear with a period Tj = = 2tt / AV deter- 
mined by the bias, with an amplitude that vanishes^ 6 - 
as 1/L. Hence, in the free case one can extract an ap- 
proximant to the steady state current fitting the Joseph- 
son oscillations with a cosine function 1 ^ of the form 
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Io(t) = I + ij cos(27rr/Tj + (p), where a = C\ or c 2 
denotes the left or right contact, and the free parameters 
of the fit are /, 1 j and (p. 

In the case of two impurities without interaction we 
find the same time scales, with minor differences. In par- 
ticular the transient time also depends on the distance 
between the two impurities, and the amplitude of the 
Josephson oscillations is also affected by R. Neverthe- 
less, as we show in Fig. |3l it is still possible to extract 
the approximant to the steady-state current by fitting 
the Josephson oscillations as mentioned above, obtaining 
an amplitude that also vanishes in the thermodynamic 
limit. In the presence of interaction, both for one and 
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seen for two different numbers of frequencies considered. 
We associate to the approximant to the steady-state cur- 
rent the uncertainty: 



AI = 



M\ 



(8) 



i=l,M 



where Tj, with i = 1,M, are the equally spaced times 
lying in the interval where the DFT is performed, /(Tj) 
is the current measured at r, and / is the reconstructed 
current. The uncertainty A J is typically within the size 
of the symbols in our plots. 

By using the procedure described above we reproduce 
in Fig. [5] the I-V characteristics of a single impurity in 
the full range of interactions and biases with excellent 
agreement with the original workS. 




Figure 4. (Color online) (a) Black continuous line: left- 
contact current I c \ for a system of L — 100 with two im- 
purities at distance R = 7, tc = 0.8io, AV = 0.5to, and 
Uc = 5to- Horizontal continuous straight line: zero frequency 
component of the DFT in the interval [20, 45] (delimited by 
vertical dashed lines), (b): DFT of the black curve in (a). 
The red dotted curve in (a) corresponds to N v = 2 frequen- 
cies: the zero frequency component and the Josephson fre- 
quency vj. The green dashed curve in (a) is found using also 
the frequencies framed by the dashed line in (b). 



Figure 5. (Color online) I-V characteristics of a system at 
tc = 0-5to with quench scheme (B). The crosses are data 
from Ref. @, the symbols are those obtained with our code 
(the parameters of our simulations are L — 96, m = 600 
states, discarded weight e < 10 -7 ). 



III. RESULTS 



two impurities, additional frequencies emerge. In Fig. [4] 
we show the current on the left contact for Uc = 5io as 
an example, where additional oscillations superimposed 
to the Josephson oscillations (they have in this case a pe- 
riod Tj ~ 12.6h/to) are clearly visible. In order to deal 
with the appearence of several frequencies, we perform a 
discrete Fourier transform (DFT) by first identifying an 
interval of time where the evolution is quasi-stationary, 
with a duration that is an integer number of Josephson 
periods Tj. Then we do a reconstruction of the current by 
picking up only the few most important frequencies from 
the DFT, which always include the zero frequency com- 
ponent (the approximant to the steady-state current), 
the Josephson frequency vj, and the one due to interac- 
tion with the highest Fourier weight v\j, as displayed in 
Fig. 21 where the quality of such a reconstruction can be 



A. Phase relations 

As is well known, the RKKY interaction is an indi- 
rect exchange interaction between two localized spins me- 
diated by the surrounding electrons of the Fermi sea£. 
In the present case, since we are dealing with spinless 
fermions, only a coupling to the density will result. The 
RKKY interaction depends on the distance between the 
impurities R via 2kp oscillations^ and is expected to in- 
duce correlations between the densities on the two dots 
and, consequently, on the currents in the contacts. We 
now show that the occupations on the dots fulfill the pre- 
dictions of the RKKY interaction. The same correlations 
are also visible in the currents, but with a phase shift. 

We take the system at half-filling, i.e. N / L = 0.5 in the 
free case and concentrate on the quasi-steady regime. In 
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the left panels of Fig. [6] we show the occupations on the 
two quantum dots. They oscillate with the Josephson 
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Figure 6. (Color online) Panels (a) and (c): time evolution of 
the number of particles on the left (blue continuous line) and 
right dot (red dashed line); panels (b) and (d): time evolution 
of the left-contact (blue continuous line) and right-contact 
(red dashed line) currents. Data for a system of L = 100, 
tc = 0.8io, AF = 0.5to, half-filling, quench scheme (B) and 
U c = 0. 

frequency uj, which characterizes also the current (see 
Sec. Ill Cl) . More interestingly we observe that when R is 
odd the densities oscillate in opposition of phase, while 
if R is even they oscillate in phase. This is a regular 
pattern which we find in all the range of R considered. 
This behavior is compatible with the 2kp oscillations of 
the RKKY interaction, as shown by Fig. [7] There it can 
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Figure 7. (Color online) Static susceptibility connecting the 
dots at sites d\ and d,2 for tc = 0.5to for half-filling. 

be seen that the static susceptibility, that displays 2fcp 
oscillations as a function of R, is positive for R odd and 
negative for R even. Therefore, for R odd the densities 
at the dots experience an effective repulsive interaction, 
while for R even it is attractive. 



If we now move to the right panels of Fig. [6] we find the 
opposite situation. When R is odd the currents oscillate 
in phase (they are exactly equal in this case) and when R 
is even they are in opposition of phase. In the latter case 
averaging the currents of the two contacts cancels out 
the oscillations. This effect is visible only in the quasi- 
stationary regime, as we can see from the left panels of 
Fig. EI 

The phase shift between densities and currents can be 
undcrtstood by noticing that when the mean density on 
a dot increases, transfer of a particle to (from) the dot 
is suppressed (enhanced) Then, for R odd, while one dot 
has a higher density, the other has a lower one. Consid- 
ering the current on the links to the left of d\ and to the 
right of d>2, charge flow is enhanced on both links when 
dl has an increased density and d2 a reduced one, while 
in the opposite case current is suppressed. On the other 
hand, when R is even, both dots have an enhanced den- 
sity or a suppressed one, such that when charge can be 
transferred on one link, the current is suppressed on the 
other. 
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Figure 8. (Color online) Same as Fig. [6] but for Uc = 5to. 



Although the evolution of the current is more involved 
in the presence of interaction due to the appearence of ad- 
ditional oscillations, the same qualitative considerations 
hold also for Uc ^ . As an example, in Fig. [8] we show 
the currents and the densities in the presence of interac- 
tion, namely at Uc = 5io . The behavior of the densities 
is very clear and analogous to the free case. However, the 
interaction enhances the amplitude of the oscillations as 
can be seen comparing Figs. [5] ancJSl In spite of the in- 
teraction, it is clearly visible that for the odd-i? case the 
currents are exactly equal and for even R an opposition 
in phase is evident. 

The phases characterizing the time evolution of the 
density and of the current described above reveal the ef- 
fect of the RKKY interaction on the slow dynamics of the 
Josephson oscillations, giving rise to sustained and con- 
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trollablc oscillations of the densities and of the currents in 
finite size systems. This fact may turn out to be observ- 
able in experiments focused on quantum dots set-ups in 
mesoscopic systems. Indeed there have been proposals 
of simulating quantum impurity systems and transport 
properties in cold atom system s 31 ' 32 . The first experi- 
mental progress done so far in this direction is the realiza- 
tion of a mesoscopic conducting channel in a cloud of Liti- 
hum atoms, performed by Brantut and collaborators^. 



B. Average current as a function of the distance 
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1. Free case 

The results shown above indicate that the dynamics of 
the current and the density is regulated by 2k-p oscilla- 
tions due to the RKKY interaction. We now investigate 
how the behavior of the steady-state current is affected 
by the distance between the impurities and the Fermi 
momentum. In Fig. [5] we show the approximant to the 
steady state current in absence of interaction as a func- 
tion of the distance between the impurities. The simplest 
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Figure 9. (Color online) Approximant to the steady state 
current as a function of the distance between noninteracting 
impurities for a system of L — 200 sites for different values of 
the contact hopping, quench scheme (B) and Uc = 0.0. Blue 
dashed line: Landauer-Biittiker formula (cf. Eq. 1111) applied 
to the case tc = 0.5£o- From top to bottom, left to right: 
AV = 0.25, 0.5, 0.8, 1.0, 1.2, 2.0t - (Fhe current is plotted in 
absolute value.). 



case is tc = to, for which the data of Fig. [9] show very 
small variations as a function of R. Indeed these varia- 
tions are only a finite size effect: as we show in Fig. [TU] 
the current for R converges in the thermodynamic 
limit to the value corresponding to R = 0. On the other 
side, from Fig. [HI we see that the curves with tc = 0.5to 
are the most sensitive to R, showing pronounced fluctua- 
tions. Furthermore, for tc = 0.5to, there is a range of R 



Figure 10. (Color online) Approximant to the steady state 
current for a system with tc = to, AV = 0.25£o and scheme 
(B) for different sizes L, as obtained from time averaging in 
the interval t = [15 -i- 41]. The dashed line is the current for 
R = 0. 



where the oscillations have the largest amplitude. This 
range changes with AV. As an example, for AV = 0.5to 
the range is given approximately by R ~ 7 -j- 19, while 
for AV = 0.75<o by R ~ 4 -f- 13. Moreover, the period of 
these oscillations is typically R — 2. 

Following the Landauer-Biittiker approac h 34 ' 35 , we 
now show that the patterns of the current of Fig. [5] can 
be understood in terms of the transmission properties 
for a single particle. Indeed, the physical mechanism at 
the root of the flow of current is that the dot, charac- 
terized by tc 7^ to, is an effective tunnel barrier with an 
energy-dependent transmission probability p s (e) (where 
the subscript s stands for single). The presence of two 
dots requires the combination of the transmission propa- 
bilities in order to compute the total probability Pd{t) 
(the subscript d standing for double). The transmission 
probability through a single dot is given byi£: 



l-e*/{Atl) 



l + ^{tl-2t c )/{At c ) 



(9) 



The total transmission probability can be obtained using 
the transfer matrix approac h 35 ' 36 and gives: 



Pd 



1 + (1 - p s ) 2 - 2(1 - p s ) cos(2fc(e)i? - 2c, 



(10) 



where <fi = kb and b is the size of the single tunnel barrier. 
In our case we have that tc is present on three sites (the 
dot and its nearest-neighbors), so b = 3. The expression 
for the combined probability cq. (|10[) is valid provided 
R > 3. Indeed, for R = 1, 2 one has to consider a single 
barrier of size 6 = 4,5 respectively. In order to obtain the 
average current, one has to integrate the transmission 
probability over the energies of current-carrying states. 
This yields^ 4 .: 



rtF+AV/2 

I(AV) = / dep d (e) . 

Je F -AV/2 



(11) 
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Our results for tc = 0.5io are the blue dashed lines of 
Fig. 13 We observe that there is a very good agreement 
with the current obtained by doing the time average. 



2. Interacting case 

We start considering the effect of a small interaction, 
namely Uc = 1-Oto- and we choose tc = 0.5to in order to 
probe if and how the interaction affects the resonances 
(Fig. [TTj) . From the comparison with the free case we 
can see first that the current is enhanced, an effect that 
becomes stronger at larger values of the bias at larger 
values of R. The enhancement of the current by interac- 
tion is also observed in the one impurity case^ (see Fig. [5] 
for small Uc < to and AV < 2io) Furthermore, it can 
be seen that the resonances observed in hte free case are 
suppressed. 




Figure 11. (Color online) Approximant to the steady-state 
current as a function of the distance between the impurities 
R for a system of L = 100 sites with Uc = l.Oto for differ- 
ent values of the bias and quench scheme (B). The currents 
in the non-interacting case (dotted, dash-dotted and dashed 
lines show the current and correspond to AV = 0.5, 0.8, 1.2to 
respectively) are shown as reference. The uncertainty on the 
value of the average current AI as discussed in Sec. Ill CI is 
within the size of the symbols in the plot. (The current is 
plotted in absolute value.) 



In Fig. [TJ] we show our results for the approximant to 
the steady-state current with increasing values of the in- 
teraction. While the current does not vary significantly 
for values of Uc lower than the band-width, a qual- 
itatively different behavior appears when Uc is larger 
than 4£o- For Uc = 6, 10to and R > 6 we interestingly 
find that the current oscillates as a function of the dis- 
tance with periodicity two, with a rather large amplitude, 
which is typical of RKKY oscillations at half-filling. The 
same behavior is also confirmed if we change the contact 
hopping, for example with tc = 0.8io- In Fig. |9] we saw 
that without interaction the current is in this case almost 
independent on R, because the single tunnel barrier has a 
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Figure 12. (Color online) Approximant to the steady-state 
current as a function of the distance between the impurities R 
for a system of L = 100 sites with tc — 0.5to and = 0.5to 
for different values of the interaction and quench scheme (B). 
In some cases the approximant to the steady-state current is 
strongly sensitive on the chosen boundaries of the stationary 
regime where the DFT is performed. The resulting different 
values of the current lie inside the error bars. (The current is 
plotted in absolute value.) 



transmission coefficient close to unity (see Eq. [9]) . Com- 
paring Figs. IT21 and [13] we find that the approximant to 
the steady-state current shows the same pattern if the 
interaction is large enough, i.e. Uc ^ 5.0to, hinting at a 
signature of the RKKY interaction. It is also remarkable 

0.6 1 — ' — i — i — i — ' — i — ' — i — i — i — ' — i — i — i 




Figure 13. (Color online) Approximant to the steady-state 
current as a function of the distance between the impurities 
R, in the presence of a small bias AV = 0.5to for a system of 
L — 100 sites for different values of the interaction and tc = 
0.8£o- The uncertainty on the value of the average current A7 
as discussed in Sec. Ill CI is within the size of the symbols in 
the plot. (The current is plotted in absolute value.) 

that the maxima of the current are of the same order 
as for the single-impurity case. This behavior and the 
features of the time evolution of the currents discussed 
in Sec. MI A"l make the large interaction regime physically 
very interesting. By putting a Fermi sea between two 



8 



interacting levels, one finds that the current can be of 
the same order of the one impurity case, but with the 
advantage of controlling the phase relations between the 
currents and the occupations. 

IV. SUMMARY 

By studying the time dependence of the current and 
the density in a one-dimensional chain in the presence of 
two interacting resonant levels, we observed the interplay 
of the RKKY interaction and the characteristics of the 
quantum dots, concerning the dynamical behavior in a 
finite system as well as the approximant of the steady- 
state current. 

Focusing on the time evolution, we found that, at fi- 
nite size, the evolution of the current in the contacts and 
the occupations of the dots are characterized by oscilla- 
tions, whose period depends on the applied bias as in the 
single dot case^, but interrelated in a way that depends 
on the size of the Fermi sea. In fact, we show that the 
densities on the dots oscillate with a relative phase which 
depends on the Fermi momentum of the Fermi sea and on 
the distance between the impurities, as expected for the 
RKKY interaction. An analogous behavior is found for 
the time evolution of the currents in the contacts, which 
are related to those of the density, but phase shifted with 



respect to them. Such a behavior could be exploited in 
experimental measurements in mesoscopic systems, giv- 
ing thus a direct access to the RKKY interaction. As 
mentioned before, experimental investigation of trans- 
port in cold atomic system o 32 ' 33 i 37 ' 38 , would be an in- 
teresting set-up, where the variations of the density in 
the quantum dots could be accessed directly. 

We have also studied the approximant to the steady 
state current, and its oscillations as a function of the 
distance between the dots. In the free case we identified 
resonances that can be traced back to the resonances 
affecting the transmission coefficients of a single particle 
propagating freely in the system. Turning interactions on 
the resonances are suppressed. However, for large values 
of the interaction we observe rather large oscillations of 
the current as a function of the distance with periodicity 
two. This corresponds to 2kp oscillations, as expected 
for the RKKY interaction 
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